Introduccion a Matlab 


Toolbox: Control System Facultad de Ingenieria-UNER 


GUIA PARA EL USO DE MATLAB PARTE 1 

GUIA DE USUARIO BASICO PARA MATLAB 
El programa Matlab 

MatLab (MATrix LABoratory) es un programa para realizar calculos numericos con 
vectores y matrices. Una de las capacidades mas atractivas es la de realizar una 
amplia variedad de graficos en dos y tres dimensiones. 

MatLab es un gran programa de calculo tecnico y cientifico, que tiene su propio 
lenguaje de programacion. Dicho lenguaje, es una herramienta de alto nivel para 
desarrollar aplicaciones tecnicas faciles de utilizar. 


Al arrancar MatLab se abre una ventana como la indicada en la Figura 1 . 


4 MATLAB 7.10.0 (R2010a) 
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La parte mas importante de la ventana inicial es la Command Windows, que aparece 
en la parte derecha (recuadrada en rojo en la Figura 2). En esta sub-ventana es 
donde se ejecutan los comandos de MatLab, a continuacion del sfmbolo caracterfstico 
(»), que indica que el programa esta preparado para recibir instrucciones. 
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Figura 2 Sub-ventanas de comandos 

En la parte izquierda de la pantalla aparecen dos ventanas tambien muy utiles: en la 
parte superior aparece la ventana (recuadrada en verde en la Figura 2) Current 
Folder, que se puede alternar con Workspace clickeando en la pestana 
correspondiente. La ventana Current Folder muestra los ficheros del directorio activo o 
actual. El Workspace contiene informacion sobre todas las variables que se hayan 
definido en esta sesion y permite ver y modificar las matrices y vectores con los que 
se este trabajando. 

En la parte inferior aparece la ventana (recuadrada en azul en la Figura 2) Command 
History que muestra los ultimos comandos ejecutados en la Command Window. Estos 
comandos se pueden volver a ejecutar haciendo doble click sobre ellos. 

La instrucciones se pueden proporcionar directamente en la venta de comando 
( Command Window) o a traves de ficheros-M (o M-files). Estos ficheros tienen la 
extension *.m y contienen conjuntos de comandos o definicion de funciones. La 
importancia de estos ficheros-M es que al teclear su nombre en la Ifnea de comandos 
y pulsar Intro, se ejecutan uno tras otro todos los comandos contenidos en dicho 
fichero. 


Aunque los ficheros *.m se pueden crear con cualquier editor de texto tal como Block 
de Notas, Notepad, Word, etc, MatLab dispone de un editor que permite crear y 
modificar estos ficheros, como ejecutarlos paso a paso para ver si contiene errores. 
La Figura 3 muestra la ventana principal del Editor. El editor muestra con diferentes 
colores los diferentes tipos o elementos constitutivos de los comandos (en verde los 
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comentarios, en rojo las cadenas de caracteres, etc.). Las Imeas de comentario se 
indican con el caracter % delante del comentario. 
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Figura 3 Ventana del editor de texto 


clc; 

close all; 
t=0: 0.1:1; 
x=2*t; 

'archivo de muestra'| 


Operaciones con vectores y matrices 

Para definir una matriz o un vector no hace falta establecer de antemano su tamano. 

Las matrices se definen por filas, los elementos de la misma fila estan separados por 
blancos o por comas, mientras que las filas estan separadas por pulsaciones intro o 
por caracteres punto y coma (;). Por ejemplo el siguiente comando define una matriz 
‘A’ de dimension 3x3: 

» A=[l 2 3; 

7 3 9] 

A = 

12 3 

■3 5 6 

7 3 9 
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De forma analoga a las matrices, es posible definir un vector fila ‘x’ en la forma 
siguiente: 

» x=[10 20 30] 

X = 


10 20 30 


En MatLab se accede a los elementos de un vector poniendo el indice entre 
parentesis, por ejemplo: 

» x ( 2 ) 
an s = 


20 

Para definir un vector columna y hay que separar los elementos por (;) o intro, 

» y= [ 11 ; 12; 13] 


V = 


11 

12 

13 


Y de la misma forma se accede a un elemento de el: 

» y{3) 


13 


MatLab tiene en cuenta la diferencia entre vectores fila y vectores columna. Si 
intentamos sumar los vectores ‘x’ e y obtendremos el siguiente mensaje de error: 

» x+y 

??? Error using ==> plus 
Matrix dimensions must agree. 

Para poder sumar los vectores necesitamos que ambos sean vectores filas o ambos 
columna. 
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Si hacemos: 

» X ’ 

an 3 = 


10 

20 

30 

Transforms el vector fila ‘x’ en un vector columns. 
Si sumamos ahora obtenemos: 

» X'lY 
an 3 = 


21 

32 

■33 

Si transformamos y en vector fila y sumamos, entonces: 

» xly 1 
ans = 


21 32 33 


Operaciones entre vectores 

MatLab puede operar con vectores por medio de operadores y por medio de 
funciones. 

Los operadores son los siguientes: 

• Suma (+) : 


» X4y' 
ans = 


21 32 33 


• Resta (-) : 

» x-y 1 


ans = 


-1 


3 


17 
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• Multiplication (*): 


1 . Vector fila por vector columna, el resultado es un numero. 

» 

an s = 


740 


2. Vector columna por vector fila, el resultado es una matriz. 


» y*x 
an s = 


110 

220 

330 

120 

240 

360 

130 

260 

390 


• Producto elemento a elemento (.*): 

1. Por un numero: 

» x . *3 

ans = 


30 60 90 


2. Por un vector elemento a elemento: 

>> >i. ‘fry 1 


an s = 


110 240 390 
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• Division elemento a elemento: 

1 . Division de todos los elementos del vector por un numero(./) 

>> x./3 

an s = 

3.3333 6.6667 10.0000 

2. Division de un numero por todos los elementos de un vector (A) 

» x . \3 

ar.3 = 

0.3000 0.1500 0.1000 

• Elevar a una potencia elemento a elemento (. A ) 

» x.“3 
ar.3 = 

1000 3000 27000 

Si anteriormente no hemos definido el vector x se hace de la misma forma, es decir: 

» [1 2 3 4]. ~2 
ans = 

1 ^ 9 16 

» [ 1234 ].*[ 1 - 11 - 1 ] 
ans = 

1-2 3-4 
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Funciones de Matlab 


Funciones Matematicas elementales: Se aplican a valores escalares o a vectores 
elemento a elemento. 


sin(x) -> funcion seno 

cos(x)-> funcion coseno 

tan(x) -> funcion tangente 

log(x) -> funcion logaritmo neperiano 

Iog10(x) -> funcion logaritmo decimal 

exp(x) -> funcion exponencial 

sqrt(x) -> funcion rafz cuadrada 

round(x) -> funcion redondeo hacia el entero mas proximo 
abs(x) -> funcion valor absoluto 


Funciones que actuan sobre vectores: Las siguientes funciones solo actuan 
sobre vectores. 


[xm,im ]=max(x) 

[ym,jm]=min(x) 

sum(x) 

mean(x) 

std(x) 


Devuelve el valor maximo xm y la posicion im del vector x 
Devuelve el valor mfnimo ym y la posicion jm 
Suma de los elementos de un vector 
Valor medio de los elementos del vector 
Desviacion estandar 


Funciones para calculos con polinomios: Para MatLab un polinomio se puede 
definir mediante un vector de coeficientes. Por ejemplo, el polinomio: 


x4-8x2+6x -10=0 


se puede representar mediante el vector [1, 0, -8, 6, -10]. MatLab puede realizar 
diversas operaciones sobre el, como por ejemplo evaluarlo para un determinado valor 
de ‘x’ 


» pol=[l 0-36 -10] 


pol = 


10-36 -10 


Si se tiene el siguiente polinomio: 


p(s) = s 3 + 3s 2 + 4 


»p=[l 3 0 4] 
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Algunas de las funciones que podemos realizar sobre este son las siguientes: 



Tambien podemos realizar: 



Graficas bi-dimensionales 

MatLab dispone de cinco funciones basicas para crear graficos 2-D. Estas funciones 
se diferencian principalmente por el tipo de escala que utilizan en los ejes de abscisas 
y de ordenadas. Estas cuatro funciones son las siguientes: 

• plot() crea un grafico a partir de vectores y/o columnas de matrices, con 
escalas lineales sobre ambos ejes. 

• loglogO idem con escala logarftmica en ambos ejes. 

• semilogx() idem con escala lineal en el eje de ordenadas y logarftmica en el 
eje de abscisas. 

• semilogyO idem con escala lineal en el eje de abscisas y logarftmica en el eje 
de ordenadas. 

• plotyyO dibuja dos funciones con dos escalas diferentes para las ordenadas, 
una a la derecha y otra a la izquierda de la figura. 
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Existen ademas otras funciones orientadas a anadir tftulos al grafico, a cada uno de 
los ejes, a dibujar una cuadncula auxiliar, a introducir texto, etc. Estas funciones son 
las siguientes: 

• title('titulo') anade un tftulo al dibujo. 

• xlabel('tal') anade una etiqueta al eje de abscisas. Con ‘xlabel off’ desaparece. 

• ylabel('cual') anade una etiqueta al eje de ordenadas. Con ‘ylabel off 
desaparece. 


Ejemplo: 


» x=[l 32353] 
x = 


13 2 3 5 3 


» plot (x) 
» 



Una segunda forma de utilizar la funcion plot() es con dos vectores como argumentos. 
En este caso los elementos del segundo vector se representan en ordenadas frente a 
los valores del primero, que se representan en abscisas. 
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Ejemplo: 


» :■:=[! 6 5 2 1 ] ; y=[l 0 4 3 1]; 
» plot (x r y} 

» 



Por ultimo, desde la ventana de graficacion pueden realizarse cambios en las formas 
de mostrar las funciones como por ejemplo modificar colores, escalas, agregar 
leyendas, etc. 


Una ventana grafica se puede dividir en m particiones horizontales y n verticales, con 
objeto de representar multiples graficos en ella. Cada una de estas subventanas tiene 
sus propios ejes, aunque otras propiedades son comunes a toda la figura. La forma 
general de este comando es: 

» siibpiot {m f n f i) 

donde m y n son el numero de subdivisiones en filas y columnas, e i es la subdivision 
que se convierte en activa. Las subdivisiones se numeran consecutivamente 
empezando por las de la primera fila, siguiendo por las de la segunda, etc. Por 
ejemplo, la siguiente secuencia de comandos genera seis graficos en la misma 
ventana: 


li 
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x=0: 1:100; 
y=sin (x) ; 
z=cos (x) ; 
w=exp (x) ; 
subplot (3, 2 , 1} 
plot £x, y) 
subplot (3 f 2 r 23 - 
plot £x, z } 
subplot (3, 2 , 3} 
plot ( x , w ) 
subplot (3,2,4)- 
plot £x, y . /x) 
subplot (3, 2 , 53- 
plot (x, sqrt (x) > 
subplot (3,2,6) 
plot (x, z . /x) 



Es posible dibujar varias funciones en una sola ventana con la instruccion: 
»hold on; 

Para desactivarla solo es necesario escribir: 

»hold off; 


t=Q : 0 . 01 : 1 
x=3*t; 
y=sin (x> ; 
p=exp (x) ; 
bold on; 
plot ( t ,x, 
plot (x, yj ; 
plot (x,p, 
bold off; 


'r') ; 

' g ' > ; 


n Figure 1 
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GUIA PARA EL USO DE MATLAB PARTE 2 
TOOLBOX DE CONTROL 

FUNCION DE TRANSFERENCIA 

MatLab es una potente herramienta para el analisis de ‘sistemas’ descriptos por 
funciones de transferencia. 

La funcion de transferencia de un sistema lineal e invariante en el tiempo, relaciona la 
transformada de Laplace de la salida con la transformada de Laplace de la entrada en 
un sistema de ecuaciones diferenciales a condiciones iniciales nulas. En forma 
generica se representa de la siguiente forma: 

Q( s \ _ 1 +—+a l s + a 0 

K U(s ) s n +b n _ 1 s n - l +...+b l s + b 0 

En sistemas reales o ffsicamente realizables m <= n. 

El polinomio del denominador igualado a cero representa la ecuacion caracterfstica 
que se utiliza ampliamente en el analisis de la estabilidad del sistema. 

Para crear funciones de transferencia en MatLab se utilizan los siguientes comandos: 


a) g=tf(num,den) 

donde “num” es un vector que contiene los coeficientes del polinomio del numerador 
de G(s) ordenado respecto a las potencias de s donde el primer elemento es el 
coeficiente que acompana a la mayor potencia de s. “den” es otro vector que contiene 
los coeficientes del polinomio del denominador de G(s) ordenados de la misma forma 
que para el numerador. 


Ejemplo de sintaxis en MatLab 
% Introducir una funcion de transferencia polinomica 


G(s) = 


s + 2.s + 3 
5 , ^+3.i' +3.^ + 1 


>> num= [ 1 f 2 f 3] ; 

>> den= [1, 3, 3, 1J ; 
>> S=tf (nuri f den) 

Transfer function: 
3~2 423+3 


s~3 + 3 s~2 +3s+l 
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b) g=zpk(z,p,k) 

Donde “z” es un vector que contienen los ceros del numerador de G(s), “p” es un 
vector que tiene los polos de G(s) y “k” es la ganancia estatica de G(s) 

Ejemplo de sintaxis en MatLab 

%Cargar en Matlab una G(s) que tiene ceros en -1 y -2, polos en -10, -3+/-3i 
% y ganancia estatica k=5 

5(5 + l)(s + 2) 

»z=[-l,-2J; (s + 10)(s A 2 + 6s + 1 8) 

» p=[-10 f -3+3i, — 3— 3i J ; 

>> h=zpk (z f p f 5) 


Zero/pole/ gain : 

5 (341) (342) 


(3+10) (s A 2 - €3 + 13) 


c) s=tf(‘s’) 


A partir de esta instruction de puede utilizar la “s” en las expresiones polinomicas de 
G(s) para que Matlab las interprete como funciones de transferencia. 


Ejemplo de sintaxis en MatLab 
% Introducir una funcion de transferencia polinomica 


G(s) = 


s + 2.s + 3 
s ^ + 3.s + 3 + 1 


» 3=tf ( 1 3 1 > 

Transfer function: 

3 

» G= ( s A 2+s*2 +3 ) / (s A 3+3*s A 2+3*s+l) 

Transfer function: 

3 a 2 - 2 3 + 3 

3^3 + 3 3^2 +3s + l. 
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Una funcion de transferencia se puede descomponer en fracciones simples, utilizando 
la instruccion ‘residue’ de MatLab. La diferencia con respecto a la forma tradicional 
de trabajo, es que, si existen polos complejos, los residuos que devuelve ‘residue’ 
corresponden a cada uno de los polos: 

7( ^) - a " ,r +a „ , - i ^ 1 ~' + --- +a ^ + a ° _ r d) | r( 2) | | r(n ) | 

U(s ) s n +b n _ 1 s n ~ 1 +...+b l s + b 0 s-p(l) s-p( 2) s-p(n ) 



Como practica ingresar la siguiente funcion de transferencia y descomponerla en 
fracciones simples : 


» n=[l 1]; 

» d=[l 3 1 7]; 

» [r,p,k]=residue(n,d) 


G(s) = 


^ + 1 

+ 3.s + s + 7 


-0.1630 

0.0815 - 0.14931 
0.0815 + 0.1493i 

P = 

-3.3307 

0.1654+ 1.4402i 
0.1654- 1.4402i 

k= [] 
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REPRESENTACION EN DIAGRAMA DE BLOQUES: 

Los diagramas de bloques representan graficamente a un sistema indicando las 
funciones realizadas por cada componente y el flujo de las senales. Generalmente los 
elementos que podemos encontrar en este tipo de representacion son: la planta 
(elemento a controlar), controladores, actuadores, sensores, entre otros. 

Un sistema sencillo de control a lazo abierto lo podemos obtener a partir de la 
interconexion en serie de la planta y del controlador. La funcion de transferencia 
obtenida para este caso en forma generica serfa: 



En ciertos sistemas pueden aparecer bloques en paralelo, y para su resolucion, 
MatLab utiliza la instruccion parallel: 



[num,den]=parallel(numl ,denl ,num2,den2) 
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En el caso que se quiera obtener la funcion de transferencia de un sistema 
realimentado podemos utilizar la funcion ‘feedback’ de MatLab para calcular la 
funcion resultante : 



Como ejemplo se calcula la funcion de transferencia total de los siguiente diagrama de 
bloques haciendo uso de las instrucciones de Matlab de dos maneras diferentes : 
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Transfer function: 

11340 s + 3240 
2 s A 2 + 1.134e004 s + 3241 


Resumen de instrucciones para la resolucion de bloques: 
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RESPUESTA TEMPORAL DE SISTEMAS: 

El objetivo que nos planteamos es el de mostrar como Matlab describe la respuesta 
temporal de cualquier sistema a diferentes entradas: 

• Entrada escalon : 


Y(t) - salida valuada en t 


num 

G(s) = — 

den 


t - tiempo al que es calculado la respuesta al escalon 


X(t) - respuesta de esatdo en t 
t - tiempo de simulacion 



[y,x,t]=step(num,den,t) 


( opcional ) 


Ejemplo usando funcion de transferencia: 


» 3=tf { ' 3 ' > 

Transfer function: 

3 

» g= 1/ (s+5) ; 

» step (g f 10) 

» 
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• Entrada impulso : 


Y(t) - salida valuada en t G(s) = t - tiempo al que es calculado la respuesta al 

den 

X(t) - respuesta de esatdo en t impulso ( opcional ) 

t - tiempo de simulacion 

[y,x,t]=impulse(num,den,t) 



Ejemplo usando funcion de transferencia: 


» g= 1/(3*345*341); 
» impulse (g r 10) 
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• Entrada arbitraria creada por el usuario : 



Ejemplo usando funcion de transferencia: 


» s=tf ( ' s 1 } ; 

» g= (s+1) / ( (3+2) *(s+3) >; 
» t=u : 0 . 1 : 100; 

» u=t *' 3 ; 

» [y f xl=l3im (g f \i r t } ; 

» plot (t r y) 

» 



21 



Introduccion a Matlab 


Toolbox: Control System Facultad de Ingenieria-UNER 


LUGAR DE RAICES: 

Teniendo un sistema como el de la 
figura la funcion transferencia a lazo 
cerrado es: 

C(s) = Gjs) 

R(s ) 1 +G(s)H(s) 

La ecuacion caracterfstica esta dada por 
igualado a cero. 



denominador de esta funcion transferencia 


G(s)H(s ) = -1 


Para que esta igualdad se cumpla, y recordando que estamos en el piano complejo, 
se deben cumplir dos condiciones: 


Condicion de angulo: 

ZG(s)H(s) = ±180°(2*+1), k = 0,1,2, 


Condicion de modulo: 


\G(s)H(s)\ = 1 

El lugar de rafces es un metodo a traves del cual se puede visualizar como varfan las 
rafces de la funcion transferencia a bucle cerrado a medida que la ganancia K del 
sistema crece de cero a infinito. Teniendo en cuenta que la ecuacion caracterfstica 
debe tener la siguiente forma: 


j + /n> + Zi)Q + z 2 )...Q + zJ 

(s+Pi)(s+p 2 ) — (s+p n ) 


Utilizando la instruccion rlocus en MatLab puede obtenerse el lugar de rafces de una 
funcion transferencia. 
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Ejemplo: 


% Siendo g la funcion en transferencia a bucle abierto. 


» g= [3+0 .5)/ (3*3* ( 3+1) +1) 

Transfer function: 

3 + 0.5 

s^3 + s A 2 + 1 

» r locus (g) 

» 



Los detalles de como trabajar con la grafica (Ifmites, tftulos, colores, etc) seran 
explicados en la seccion de respuesta en frecuencia ya que son similares a los que se 
pueden realizar sobre los diagramas de Bode. 

RESPUESTA EN FRECUENCIA: 

La funcion respuesta en frecuencia representa la manera en que el sistema responde 
a entradas de diferentes frecuencias. Para calcularla, se reemplaza la variable ‘s’ de la 
funcion de transferencia por ‘jco’. 

Bode 

La orden ‘bode’ calcula las magnitudes y los angulos de fase de la respuesta en 
frecuencia de sistemas continuos, lineales e invariantes en el tiempo. Los diagramas 
de Bode se utilizan frecuentemente para analizar y disenar sistemas de control. Estos 
diagramas indican el margen de ganancia, el margen de fase, la ganancia, el ancho 
de banda, etc. 
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Cuando se introduce la orden ‘bode’ (sin argumentos del lado izquierdo), MatLab 
realiza el diagrama de bode sobre la pantalla, permitiendo ademas recuperar datos de 
los diferentes puntos. 

El diagrama de bode, resulta de mucha utilidad a la hora de establecer la estabilidad 
relativa de los sistemas. Para esto se utilizan los margenes de estabilidad: 

Margen de Ganancia: Es valor de la ganancia en db cuando la fase cruza por -180° 
cambiado de signo. Para la estabilidad, este parametro debe ser mayor a 6 db. 

Margen de Fase: Es el valor que toma la fase cuando la Magnitud cruza por Odb 
sumado a 180. Para la estabilidad, este parametro debe ser mayor a 30°. 


» g 

Transfer function: 
3 + 0.5 


s^3 + 3^2 + 1 

» bode [g) 

» 



En la grafica se muestran los valores correspondientes al margen de fase y al margen 
de ganancia del sistema. Para que estos indicadores aparezcan en pantalla, se debe 
hacer click con el boton derecho del mouse sobre la grafica, seleccionar la opcion 
‘Caracterfsticas’ y luego elegir ‘Margenes de estabilidad rmnimos’. Hacer click con el 
boton izquierdo del mouse sobre los puntos para que aparezca el letrero. 

Al hacer click con el boton derecho del mouse sobre las graficas aparecen varias 
opciones: 

• System: indica con que color se estan representando cada funcion de 
transferencia en caso de que estemos graficando varias de ellas. De igual 
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modo, se puede seleccionar para que solo aparezca en la grafica la funcion de 
transferencia que a rrn me interesa ver. 

• Caracteristicas: Posee tres opciones, de las cuales los dos ultimas tiene 
relacion a los margenes de ganacia y fase. La primer opcion, permite encontrar 
un pico en la magnitud en caso de que la grafica lo posea. 

• Show: Permite seleccionar si se quiere visualizar la magnitud, la fase o ambas 
graficas en la misma ventana. 

• Grid: Agrega una grilla que facilita la observacion de las graficas. 

• Full view: Optimiza la visualizacion de las graficas a fin de utilizar todo el 
espacio disponible con el trazado de la grafica. 

• Propiedades: 



Dicha ventana permite modificar varios parametros de la grafica como por ejemplo: 
etiquetas sobre los ejes, Ifmites de la grafica, unidades de las variables, estilos de 
grafica ( semilogarftmico, lineales) y desde esta ventana tambien se pueden varias las 
opciones antes descriptas. 

Para obtener informacion de la grafica o de cualquier punto de la misma, se debe 
hacer click con el boton izquierdo sobre la grafica. A continuacion aparecera un punto 
cuadrado negro, con la informacion de dicho punto. Si colocamos el mouse sobre el 
punto y mantenemos apretado el boton izquierdo, podemos arrastrar el punto sobre la 
grafica hasta cualquier punto de interes de la misma. 

Al escribir la sentencia ‘margin (funcion de transferencia)’ se obtiene lo siguiente: 
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Por otra parte, es posible obtener los margenes de ganancia y los margenes de fase 
mediante la funcion ‘margin’. 

[mg,mf,wmg,wmf]=margin(num, den); 

Mg es el margen de ganancia del sistema. 

Mf es el margen de fase. 

Wmg es la frecuencia a la cual sucede el cruce de la fase por -180 Q . 

Wmf es la frecuencia de cruce de ganancia por cero. 


Otra forma de introducir la sentencia bode en Matlab es la siguiente: 

[mag, fase, w] = bode(num,den,w) 

Las matrices mag y fase contienen las magnitudes y los angulos de fase de la 
respuesta en frecuencia del sistema evaluados en los puntos de frecuencia 
especificados por el usuario. 

El vector w se define, por ejemplo, de la siguiente manera: 

w = logspace(-2,3,100) 

Siendo el valor representado por -2, en este caso, la potencia del Ifmite menor de la 
grafica. El 3 representa la potencia del Ifmite superior. Por ultimo, el 100 es la cantidad 
de espacios separados logarftmicamente entre los extremos. 
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Nyquist 

Es otra manera de representar la respuesta en frecuencia del sistema. Consiste en 
una grafica donde el eje ‘x’ se denomina eje real y el ‘y’, eje imaginario. La imagen se 
construye con la magnitud y la fase segun como varien con la frecuencia u>. 


» g=l/ (s+l) ; 
» nyquist (g) 
» 



Por otra parte, es posible obtener los margenes de ganancia y los margenes de fase 
mediante la funcion ‘margin’. 

[mg,mf,wmg,wmf]=margin(num, den); 

Mg es el margen de ganancia del sistema. 

Mf es el margen de fase. 

Wmg es la frecuencia a la cual sucede el cruce de la fase por -1 80 Q . 

Wmf es la frecuencia de cruce de ganancia por cero. 
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